NVU perspective on simple liquids' quasiuniversality 
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The last half century of research into the structure, dynamics, and thermodynamics of simple 
liquids has revealed a number of approximate universalities. This paper argues that simple liquids' 
reduced-coordinate constant-potential-energy hypersurfaces constitute a quasiuniversal family of 
compact Riemannian manifolds, which is parameterized by a single number. From this follows the 
previously discussed quasiuniversalities. 



I. INTRODUCTION 

By their contrast to complex liquids simple liq- 

uids are traditionally defined as systems of spherically 
symmetric particles interacting via pair forces [6l-[l6j. A 
number of such systems, however, like the Gaussian core 
model (T7l - [20j . the Lennard- Jones Gaussian model 2lj], 
the Jagla and related discontinuous-force models [22H25j , 
and others [25l - [29j exhibit quite complex behavior. On 
the other hand, van der Waals molecular liquids are gen- 
erally regular and "simple" in their properti es (H , [Til ] . In 
view of these facts we recently suggested [3Q] defining 
instead liquid simplicity from the property of strong cor- 
relations between equilibrium virial and potential energy 
fluctuations in the NVT ensemble [3ll - l35| . This is how 
the term "simple liquid" is used below. In practice there 
is considerable overlap between the two definitions, for 
instance the Lennard- Jones liquid and related systems 
are simple in both senses. One notable difference, how- 
ever, is that realistic liquids are only simple in the present 
meaning of the term in part of their phase diagram - 
simplicity does not apply near or at the critical point or 
at supercritical and gas states (where different kinds of 
simplicity apply, of course). With regard to real liquids, 
it appears that most or all van der Waals bonded and 
metallic liquids are simple, whereas covalently bonded, 
hydrogen-bonded, and strongly ionic liquids are gener- 
ally not simple because competing interactions weaken 
the virial potential-energy correlations |30| . 

Simple liquids are characterized by having isomorphs 
in their phase diagram 34]. An isomorph is an equiva- 
lence class of the following equivalence relation: two state 
points are isomorphic if all pairs of physically relevant 
microconfigurations of the state points, which trivially 
scale into one another, have the same configuration-space 
canonical probability. Only inverse-power-law (IPL) po- 
tentials have exact isomorphs, but all strongly correlating 
liquids have isomorphs to a good approximation [34j . 

Simple liquids' simple properties derive from the fact 
that the existence of isomorphs implies that the phase 
diagram is effectively one-dimensional instead of two- 
dimensional for the several properties that are isomorph 



invariant. Examples of such properties are [34|: Newto- 
nian and Brownian reduced-unit dynamics, reduced-unit 
static structure factors of any order, the excess entropy, 
the isochoric heat capacity. Another example is that 
for any simple liquid melting defines an isomorph in the 
phase diagram. This implies invariance along the melting 
curve of, e.g., excess entropy, reduced viscosity, reduced 
heat conductivity, reduced diffusion constant, etc, as well 
as invariance of the Lindemann melting criterion |3~I1 . |35| . 

The isomorph theory explains a number of previo usly 
noted regularities relating to any given simple liquid [34| . 
The theory, however, cannot explain the several intrigu- 
ing similarities between different simple liquids known for 
a long time. This is the focus of the present paper, which 
views simple-liquid quasiuniversality from an NVU per- 
spective. 

NVU dynamics [36j |. which is inspired by earlier ap- 
proaches to dynamics conserving [37l |38| or limiting 
[39( the potential energy, is defined as geodesic mo- 
tion on the constant-potential-energy hypersurface. If 
R = (ri, rjv) is the 3A-dimensional position vec- 
tor describing a system of N particles, R^ the posi- 
tion vector at time step i, and the corresponding 
3iV-dimensional force vector, the NVU algorithm [3(| is 
R s:+ i = 2R i -R i _ 1 -2[F i -_(Ri - Ri_i)]Fi/F?. As shown 
in Ref. |4fj, if m is the particle mass and At the time step 
of the Verlet algorithm R. i+ i = 2Rj -Rj_i +F l (At) 2 /m, 
the NVU and Verlet algorithms are equivalent in the 
thermodynamic limit because the fluctuations of the 
NVU- force prefactor become insignificant as N — s- oo. 
Consequently, the radial distribution function, diffusion 
constant, coherent and incoherent intermediate scatter- 
ing functions, etc, are identical in the thermodynamic 
limits of NVU and standard NVE or NVT Newtonian 
dynamics. This has been confirmed in computer simula- 
tions of atomic as well as molecular models [3(1 H3, 5H ■ 

NVU dynamics provides an alternative view of a liq- 
uid's molecular dynamics. At any given state point all 
information about the liquid's structure and dynamics 
is encoded in its constant-potential-energy hypersurface 
O. If (U) is the average potential energy at the state 
point, this compact Riemannian differentiable manifold 
is defined by 



* dyre@ruc.dk 



n = {R\U(R) = {U)}. (1) 
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il is a so-called level set of the function U(R). If stan- 
dard periodic boundary conditions are employed, ft is 
embedded as a (37V— l)-dimensional hypersurface in the 
3iV-dimensional torus. The manifold ft is not only de- 
fined for simple liquids, of course, but for all liquids - in 
fact also for all solids and gasses. 

The present paper considers systems of N identical 
particles in volume V interacting via pairwise additive 
forces. Thermodynamic quantities are excess quantities, 
i.e., in excess of the corresponding ideal gas quantities at 
the same density and temperature. Thus S is the excess 
extensive entropy (S < 0) and Cy the excess extensive 
isochoric specific heat, which we refer to as just "entropy" 
and "specific heat" . The corresponding intensive quanti- 
ties are denoted by lower-case letters, i.e., s = S/N and 
cy = Cy/N. Reduced quantities are marked by a tilde. 

It is clear from the results of many years of research 
into the liquid state that there is no such thing as abso- 
lute universality, even among narrowly restricted classes 
of liquids. Thus any theory predicting liquid-state uni- 
versality invariably breaks down because it is too simple. 
The philosophy of this paper is that quasiuniversalities 
provide useful insights, even though these are known from 
the outset to be only approximate. In this connection it 
is an obvious conjecture that genuine simple-liquid uni- 
versality is approached as the spatial dimension increases 
towards infinity - if this is the case, we have for liquids 
a situation reminiscent of that of critical phenomena. 

Section [TT] argues for the existence of a quasiuniver- 
sal family of constant-potential-energy hypersurfaces for 
simple liquids, parameterized by just one parameter. The 
argument presented is not rigorous, but meant to sug- 
gest one possible route for justifying quasiuniversality. 
The reader may choose to accept quasiuniversality in its 
NVU formulation as a postulate, skip Sec. lU and pro- 
ceed to Sec. IIIII that derives and discusses a number of 
consequences of the NVU formulation of simple liquids' 
quasiuniversality. Section IIVI returns to the question of 
what causes quasiuniversality. Finally, Sec. |V] gives a 
brief summary. 



II. A SINGLE-PARAMETER FAMILY OF 
HYPERSURFACES COMMON TO ALL SIMPLE 
LIQUIDS 

Any state point in the phase diagram of a liquid gives 
rise to a constant-potential-energy hypersurface fi as de- 
fined in Eq. (TT]). At first sight these manifolds may ap- 
pear to be completely characterized by the single number 
([/), but actually the system volume V is a second pa- 
rameter that is implicit in the definition of U (R) . Thus 
for liquids in general, Q is described by two parameters, 
corresponding to the two dimensions of phase space. An- 
other way of seeing this is to write the pair potential as 
v (r) = e<j)(r/a); it is then clear that any given state point 
gives rise to two dimensionless numbers characterizing 
the state point's reduced-coordinate constant-potential- 



energy hypersurface fl, namely <r 3 p (where p — N/V is 
the particle density) and e/UbT. 

Using reduced units means meas uring le ngth in units 
of p -1 / 3 , time in units of p -1 / 3 \j 'm/hsT , and energy 
in units of fc^T. The reduced 3iV-dimensional position 
vector is thus defined by R = p 1 / 3 R. Appendix A of 
Ref. showed that a liquid has strong virial potential- 
energy correlations (i.e., is simple) if and only if the liquid 
has good isomorphs, and that this happens if and only 
if the liquid's reduced-unit constant-potential-energy hy- 
persurfaces are invariant along certain curves in the phase 
diagram (the isomorphs). Thus for any given simple 
liquid a single number, A, parameterizes the reduced- 
coordinate constant-potential-energy hypersurfaces. We 
indicate this by writing 

O = f2(A) (2) 

where 

O = {R| Uip-^R) = (U)} . (3) 

The isomorph theory says nothing about how these hy- 
persurfaces compare between different simple liquids. We 
now argue that the family f2(A) is quasiuniversal, i.e., ap- 
proximately the same for all simple liquids. 

The only systems with 100% correlation between NVT 
equilibrium fluctuations of virial (W = — 1/3 J2i r i ' Vji7) 
and potential energy are the inverse-power-law (IPL) sys- 
tems, for which the potential energy scales with inter- 
particle distance as oc r~ n . For these systems the 100% 
correlation follows from the identity U(XR) = A~™C/(R) 
and Euler's theorem for homogeneous functions. A sim- 
ple liquid by definition has strong WU correlations. This 
implies that at each of its strongly correlating state points 
the liquid behaves much like an IPL system. The value of 
the effective IPL exponent n generally varies with state 
point. For Lennard- Jones (LJ) liquids at typical state 
points n is fairly constant, between 5 and 6. For the 
99.9% correlating "repulsive LJ liquid" defined by the 
pair potential e[(r/cr)~ 12 + (r/cr)~ 6 ] [42j, the exponent 
varies from n = 12 at high densities to n = 6 at low 
densities. How much the effective IPL exponent varies 
throughout the phase diagram, however, is of no impor- 
tance for the below arguments for 17(A) quasiuniversality. 

IPL systems with different exponents have very simi- 
lar ("quasiuniversal") behavior. This was noted a long 
time ago in relation to these systems' structure and DC 
dynamic properties expressed, e.g., via the diffusion con- 
stant [43T - [51j . During the last decade IPL quasiuniversal- 
ity has come into focus again [52-54J and been extended 
to include general dynamic properties, termed "dynamic 
equivalence" by Medina- Noyola and coworkers. Dy namic 
equivalence has been established for Brownian [55|, [56| as 
well as Newtonian [5TM60| dynamics. 

Because structure and dynamics are both encoded in 
il, IPL quasiuniversality follows if all IPL systems have 
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almost identical constant-potential-energy hypersurfaces. 
But why should this be the case? To understand this, 
we consider two infinitesimally close configurations with 
same n-IPL potential energy and show that they for all m 
to a good approximation have the same m-IPL potential 
energy. 

The m-IPL potential energy is given by U m = 
£mJ2ij( r ij/ a )~ m ' m w bich Tij is the distance between 
particles i and j. The change in U m between two nearby 
configurations is given by 



5U m = —me, 



p)"""' Snj 



(4) 



From this it follows that 
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The factor ln^ry/* 7 ) in Eq. ([5]), which does not vary 
much because it is a logarithm, is multiplied by the fac- 
tor (r^/cr) - " 1-1 that varies a lot. The product will be 
dominated by a rather narrow range of interparticle dis- 
tances, the most important of which is denoted by (r) m , 
a quantity that depends on state point and m. Asa good 
approximation one can replace ln(ry/ <7 ) by ln( (r) m /o~), 
and Eq. now becomes 



_d_ 

dm 
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me r> 
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me r 



(6) 



Recall that SU n = 0. Since SU m /(me m ) — is the unique 
solution to this first-order differential equation that obeys 
SU n = 0, it follows that SU m — for all m. - Note 
that the approximation made in replacing ln(ry-/(r) by 
a constant is questionable when m < 2, in which case 
no narrow range of interparticle distances dominates Eq. 
(O because J2j r ij™~ 1 diverges. 

Within the above approximation the potential-energy 
functions U n (R) have the same constant-potential- 
energy hypersurfaces for all n > 2. Possibly, this applies 
for all n. Since a simple liquid at a given state point 
may be approximated by an IPL system, this means that 
all simple liquids have approximately the same constant- 
potential-energy hypersurfaces. In other words, to a good 
approximation a single family of manifolds parameterized 
by one parameter, £1(A), is common to all simple liquids. 
Each isomorph of a simple liquid corresponds to a par- 
ticular value of A, i.e., to one specific manifold O. 



III. CONSEQUENCES OF fi(A) 
QUASIUNIVERSALITY 

Different IPL systems exhibit close similarities with re- 
spect to structure and dynamics, similarities that extend 



to all other simple liquids Because reduced-unit 

structure and dynamics are both encoded in O, any two 
state points of two different IPL systems with same O 
have the same structure and dynamics. The extension of 
similarities to all simple liquids follows from 0(A) quasi- 
universality. Note that the hard-sphere (HS) system also 
exhibits these quasiuniversalities because it is the n — > oo 
limit of the ri-IPL system. 

The Young- Andersen approximate scaling principle 
|U [Hy. This principle states that if two liquids at two 
state points have the same reduced-unit radial distribu- 
tion function g(r), they have the same reduced- unit dy- 
namics. g(f) is determined from O, so having the same 
g(f) implies having the same O. This implies same dy- 
namics. 

Quasiuniversality of the order-parameter maps of 
Debenedetti and coworkers ff^iVldj. Plotting a trans- 
lational order parameter versus an orientational or- 
der parameter for various state points leads to a one- 
dimensional curve for any simple liquid, because both 
order parameters are isomorph invariant [34 |. The ap- 
proximate identity between order-parameter curves of 
different simple liquids follows from 0(A) quasiuniver- 
sality, since O determines both order parameters. 

Excess entropy scaling ffJ^T^]. Rosenfeld noted in 
1977 that the reduced-unit diffusion constants D of dif- 
ferent simple liquids have an approximately universal de- 
pendence on the entropy per particle, s [74| . This type of 
quasiuniversality applies also, e.g., for the heat conduc- 
tivity as a function of excess entropy [78[ . For any simple 
liquid, since both D and s are isomorph invariant, one 
quantity is a function of the other. Quasiuniversality of 
the function D{s) is a consequence of the fact that D 
and s are both encoded in O (the entropy S — Ns is the 
logarithm of the area of O). 

The Lindemann melting criterion According 
to the Lindemann criterion a crystal melts when the vi- 
brational mean-square displacement obeys yj (x 2 ) ~ 0.1 
where x — xp 1 ^ 3 in which x is the atomic vibrational dis- 
placement from equilibrium in an axis direction. The 
melt ing curve in the phase diagram is an isomorph 
[HI l35| . so melting takes place for a particular mani- 
fold O c of the crystalline state. This manifold determines 
(x 2 ). Thus any simple crystal melts when (x 2 ) reaches a 
certain, quasiuniversal value. The Lindemann criterion 
and generalizations [7^, [H| have been criticized on 
the grounds that they are single-phase criteria, whereas 
melting occurs when the crystal and liquid free energies 
are the same, so any melting criterion should involve both 
phases. One possible resolution of this paradox is that 
the Lindemann criterion does not, in fact, determine the 
melting line, but a spinodal at a slightly higher temper- 
ature where the crystal becomes mechanically unstable 
[H, |83|. Alternatively, for the class of simple liquids, 
O quasiuniversality implies that there is basically just 
one melting process, which takes place at the state point 
where the specific crystalline manifold is O c . Any single- 
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phase melting criterion referring to this manifold applies 
universally for simple liquids. 

Melting rules referring to the liquid. Quasiuniversality 
of these follow from the fact that SI is quasiuniversal also 
on the liquid side of melting. For instance, this implies 
the Hansen- Verlet rule that a liquid crystallizes when the 
first peak of the radial distribution function reaches the 
value 2.85 [44(. Likewise, any simple liquid's cy is close to 
3ks at freezing [8^, |86[ (it is shown below that cy is en- 
coded in SI). Other quasiuniversal melting rules similarly 
follow from S1(A) quasiuniversality |82j. Examples are the 
Andrade equation from 1934 predicting a quasiuniversal 
value of the reduced- unit melting point viscosity [87l l88j , 
the Raveche-Mountain-Street criterion (89j of a quasiuni- 
versal value of the ratio between the values of the maxi- 
mum and minimum of the radial distribution function at 
freezing, Lyapunov-exponent based criteria [83| . or the 
criterion of zero higher-than-second-order liquid config- 
urational entropy at crystallization [90|. Note that the 
theory predicts a quasiuniversal constant- volume melting 
entropy for simple liquids, which is consistent with exper- 
iment |85l. l9~Ll| . - We finally mention the Egelstaff-Widom 
quasiuniversality [g^] , according to which the isothermal 
compressibility times the liquid surface tension is almost 
constant for liquids close to the triple point, equal to a 
few tenths of Angstroms. By the compressibility theorem 
the reduced isothermal compressibility is an isomorph in- 
variant. Thus the Egelstaff-Widom quasiuniversality fol- 
lows if the reduced surface tension is also an isomorph 
invariant, which is likely the case. 

Algebraic closedness properties of the class of simple 
potentials. Suppose Ui(R) and U 2 (R) are both poten- 
tials of simple liquids, i.e., strongly correlating. Then 
their sum t/i(R) + U 2 (R) and product f7 1 (R)[/ 2 (R) are 
also strongly correlating potentials. This does not fol- 
low merely from the property of strong correlation be- 
tween virial and potential energy NVT equilibrium fluc- 
tuations, but from the quasiuniversality of simple liquids: 
Since Ui(R) and U 2 (R) are both constant on the man- 
ifolds f2( A) , this applies also for their sum and product. 
In particular, note the following property. Writing a sim- 
ple pair potentials as v(r) = e4>(r/a), the derivative with 
respect to er, dv(r)/dcr, is also simple. 

Additivity of thermodynamic quantities. Suppose 
U(R) = t/i(R) + U 2 (R) in which U t (R) and U 2 (R) de- 
fine a simple liquid; an example is when U (R) is the LJ 
potential, Ui(R) is an n = 12 IPL function, and U 2 (R) is 
minus ann = 6 IPL. Then U(R) defines a simple liquid, 
and as functions of density and entropy the correspond- 
ing temperatures obey T(p, S) = T x (p, S) + T 2 {p, S) [I| . 
To show this, note that since the entropy determines SI, 
at given values of p and S the three constant-potential- 
energy manifolds are identical: Q = fij = fl 2 . This 
implies that U(p,S) — U\{p, S) + U 2 (p,S), from which 
the required T{p,S) = T^p, S) + T 2 {p,S) follows from 
the definition of temperature T — (dU/dS) p . The ther- 
modynamic relation W — (dU/dp)s similarly implies 
additivity of virials: W(p,S) = Wx (p, S) + W 2 (p, S) . 



The (excess) Helmholtz free energy F, (excess) Gibbs 
free energy G, and (excess) enthalpy H are likewise 
also additive: F(p,S) = F^S) + F 2 (p,S), G(p,S) = 
G 1 (p,S) + G 2 (p,S), H(p, S) = H 1 (p,S) + H 2 (p,S). 

An application is "additivity of melting temperatures" 
[49j : Since crystallization for all simple liquids takes place 
at a certain value of the liquid entropy, at any given den- 
sity one has T m = T m _\ + T m ^ 2 . An IPL liquid's melting 
temperature scales with density as T m oc p n ' 3 , so for the 
LJ liquid an expression of the form T m — Ap 4 — Bp 2 
applies |9|. 

A quasiuniversal equation of state. It was recently 
shown that simple liquids have simple thermodynamics 
in the sense that temperature factorizes into a prod- 
uct of a function of entropy and a function of density, 
T = f(s)h(p) 0. We now show that the function /(s) 
is quasiuniversal. This justifies writing the equation of 
state as 

T = f (s)h(p) . (7) 

The point is that the specific heat, like the entropy, de- 
pends only on the geometry of S7 and not, for example, 
on the specific details of the potential or the tempera- 
ture. This can be shown by first writing cy in terms of 
fluctuations of canonical ensemble probabilities, and then 
relating the latter to microcanonical (NVU) probabilities. 
According to Einstein C v = ((AC/) 2 ) /k B T 2 in which the 
average refers to the canonical ensemble. In terms of 
the configuration-space probability p oc exp(— U/ksT) 
this may be written Cy = fcs((Alnp) 2 ). The canonical 
ensemble is realized from the microcanonical NVU en- 
semble in the usual textbook way by considering a small 
subvolume V m of the total volume V. On average V m 
contains m particles where m/N = V m /V. Each config- 
uration of m particles in V m , (rj., ...,r m ), has a proba- 
bility p(ri, r m ) that can be calculated from the man- 
ifold SI (or SI) by integrating out the remaining degrees 
of freedom. In the microcanonical ensemble this corre- 
sponds to the volume of the lower-dimensional cut of S7 
corresponding to fixing the positions of particles 1 to m, 
and p(ri, r m ) is therefore a purely geometric quan- 
tity. More accurately, we argue as follows. The set of 
configurations in SI with precisely m particles in volume 
V m is denoted by Sl m c SI. Integrating out the remain- 
ing degrees of freedom from the configurations in Sl m 
determines p(ri, r m ), so this function is given by S7 
(or SI). From cy = fcs((Alnp) 2 )/m it now follows that 
cy = cy(X). Since cy ~ (ds/d\nT)p and s = s(A), 
this implies that at fixed density d\nT — cj)(X)dX for 
some function </>(A). This means that, while for two sim- 
ple liquids the temperatures corresponding to the same 
manifold SI may differ, the relative temperature changes 
at fixed density between different Sis are the same. By 
integration this implies that for each simple liquid one 
can write T = $(A)T»( ( o). Combining this with the 
equation of state T — f(s)h(p) shows that the function 
f(s) is determined by A, i.e., by the manifold SI. Thus 
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f2(A) quasiuniversality implies that f(s) is quasiuniver- 
sal: f(s) — fo{s). The function T*(p) — h{p) is not 
quasiuniversal; it reflects how the liquid's characteristic 
energy scale depends on density [95[. 

Quasiuniversality of simple liquids' specific-heat tem- 
perature dependence. Eliminating A between cy(A) and 
T = fo(s(X))h(p) leads to c v = F Q (T/h{p)) for a qua- 
siuniversal function F . This is consistent with the 
Rosenfeld-Tarazona expression cy oc T -2 / 5 [75[, which as 
shown by computer simulations applies to a good approx- 
imation not only for all IPL system s, bu t also for L J-type 
liquids and other simple liquids |96Ml0l| . Note that cy cx 
T~ 2 / 5 implies s cx -T~ 2 / 5 since c v = (ds/d\nT) p and 
s -> for T -> oo. This means that / (s) oc (-s) -5 / 2 . 

Quasiuniversal isochoric fragility of simple liquids. 
The reduced-unit relaxation time f is determined by f2. 
This not only means that f is a unique function of the 
excess entropy ("excess entropy scaling"), it also im- 
plies a quasiuniversal temperature dependence of f at 
constant density: The quasiuniversal equation of state 
Eq. ([7]) implies that at any given density, entropy is a 
quasiuniversal function of temperature: s = s(T/h(p)). 
This implies quasiuniversality of the form f = f(T/h(p)) 
[9o| . In particular, at co nstan t density Angell's fragility, 
— dlog 1Q (f)/dhiT\T=T g jl02l |. is a quasiuniversal num- 
ber for a given cooling ra te defining the glass transition 
temperature T g [To34l04 |. This is reminiscent of the uni- 
versal temperature depend ence of the viscosity discussed 
in 1996 by Kivelson et al. [105| . although these authors 
subtracted the high-temperature activation energy before 
demonstrating data collapsed. A quasiuniversal isochoric 
fragility is consistent with simulations of De Michele et 
at, who foun d that different IPL systems have the same 
fragility |l06l |. The prediction is not entirely consistent 
with experiment, although there does seem to be a ten- 
dency that van der Waals liquids have isochoric fragilities 
not far from 50 [107| | ; in this connection it is worth notic- 
ing that it is an experimental challenge to determine the 
isochoric fragility accurately. 

The hard-sphere system. Temperature plays no role for 
the configurational degrees of freedom of the HS liquid. 
Since only density is important, the HS phase diagram 
is effectively one-dimensional. This brings to mind iso- 
morphs, the existence of which implies that a simple liq- 
uid's phase diagram is also effectively one-dimensional. 
Is the HS liquid simple? Since its potential energy is 
identically zero whereas the virial is not, the HS liquid 
is not simple in the sense of the term used here. In our 
opinion, the HS liquid should be thought of more as the 
n — > oo limit of an n-IPL system than as a physical liquid 
of its own right. When a simple liquid is modeled by a HS 
system, each of the liquid's isomorphs corresponds to one 
value of the HS packing fraction n. This establishes a one- 
to-one correspondence X <-> rj, which explains why simple 
liquids' entropy, relaxation time, viscosity, etc, have been 
found to be quasiuniversal functions of the n parameter 
of the HS "reference" system. 

Arguments for quasiuniversality were traditionally 



based on the fact that any simple liquid is well rep- 
resented by the HS reference system @, [H, [Toi - tTTa ] . 
In the standard view, the HS system is useful because 
it captures the essence of liquids' harsh repulsive forces 
[llOl lllll Ill3j | . This picture is appealing, but runs into 
problems with known facts. On the one hand, IPL qua- 
siuniversality extends down to n = 3 o r n = 4, in fact 
for some quantities down to n = 1 [74], Ill4| , where the 
repulsions are quite smooth. On the other hand, there 
are several systems with harsh repulsive forces that ex- 
hibit anomalous behavior not captured by the HS system 
[HHHI. From the NVU perspective, the HS system's 
usefulness is not the explanation of simple liquids' quasi- 
universality, but a consequence of quasiuniversality: since 
all n-IPL systems are simple and quasiuniversal, the HS 
system inherits their simple properties because it is the 
n — > oo limit of an n-IPL system. 

Entropy's role. Theories relating entropy to a liquid's 
relaxation time go back in time at least to Bestul and 
Chang, who in 1964 noted that the glass transitions of 
different glass-forming liquid s occ ur at virtually the same 
value of the (excess) entropy [1151 ] . Since the glass transi- 
tion for a given cooling rate takes place when the liquid's 
relaxation time reaches a certain value, by generalization 
to other cooling rates this result implies that the relax- 
ation time is a (quasi) universal function of entropy. Inde- 
pendently, based on computer simulations and analytical 
arguments, Rosenfeld in 1977 proposed excess entropy 
scaling [74|. These two results, as well as the Adam- 
Gibbs mo del from 1965 in which entropy is also crucial 
may appear counterintuitive since entropy is 



global property: How can a global property control the 
relaxation time, which is determined as an average of 
local properties? For simple liquids il(A) quasiuniversal- 
ity provides a straightforward answer. Entropy identifies 
the relevant manifold f2, which determines the relaxation 
time. Accordingly, other unique markers of should be 
equally useful for "determining" the relaxation time, for 
instance the two-particle entropy that Dzu guto v in 1996 
suggested controls the relaxation time [lU, Ill9l ] . 

Characterizing £1 via the mean curvature. It is difficult 
to visualize a high-dimensional differentiable manifold. 
A primitive analog is a two-dimensional closed surface 
in ordinary three-dimensional space. Such a surface has 
two obvious characteristics, its area and its mean curva- 
ture. The latter is conveniently quantified in terms of the 
so-called (average) radius of curvature 91, defined as the 
inverse of the curvature. Both the area and the curvature 
conce pts gener alize to multidimensional Riemannian sur- 
faces [l20Ml22j | ; fortunately pedagogical introductions to 
this fasci nating branch of mathematics are now available 
HM23. The entropy is the logarithm of the mani- 



fold's area, but what is the physical interpretation of the 
radius of curvature? To answer this we start fro m th e 
confi guration-space canonical ensemble expression [l26l — 

D3, 
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knT 



((W) 2 ) 
(V 2 l7) 



(8) 



which is derived by partial integration of (V 2 £/) oc 
/ dRV ■ Vf/(R)exp[-£/(R)/fc B T]. Because of ensem- 
ble equivalence, the configuration-space microcanonical 
ensemble may be used to calculate the averages in Eq. 

as integrals over fL The inverse radius of curvature 
at a poin t on a d-dimensional hypersurface f2 is V • n/d 
[I24lll25| where n is the normal vector to f2 at the point, 
n = VU/|VU|. To leading order in l/VN variations in 
the denominator are insignificant, so this implies for the 
average inverse radius of curvature 1 /£H (replacing 3 AT— 1 
by3iV) 



1 (V 2 LQ 
3N (|VC/|) 



(9) 



If the average length of the 3^- dimensional force vector 
F = — VU is denoted by F, because fluctuations are in- 
significant as N — > oo, one has F 2 = ((VJ7) 2 ) and Eqs. 
© and © imply 



fRF = 3Nk B T. 



(10) 



A small radius of curvature of f2 thus corresponds phys- 
ically to a large average force. As N — > oo, (V 2 C/) ~ N 
and (|V?7|) ~ \fN which implies D\ ~ y/N . Likewise, 
F ~ \/N as A — > oo. If one defines <K = p^fR/VSN 



and F 



-1/3 



these two dimensionless 



quantities are independent of N in the thermodynamic 
limit and related by 



1 . 



(11) 



Thus the manifold f2's curvature is basically F, a quantity 
that provides an alternative to the entropy for character- 
izing f2. 

The entropy and F both have simple geometric inter- 
pretations, but the curvature F has the advantage of be- 
ing the average of a locally defined quantity. A further 
advantage is that, since fluctuations are unimportant in 
the thermodynamic limit, F may be calculated from a 
short-time simulation. Moreover, F may be calculated 
from standard NVE or NVT simulations. We conclude 
that the quasiuniversal entropy dependencies observed 
for simple liquids' structure and dynamics may equally 
well be interpreted as quasiuniversal dependencies on F, 
the curvature of fL Interestingly, this is the quantity that 
controls the relaxation time in the entropi c barrier hop- 
ping theory of Schweizer and co-workers [l29Nl31 1 (ex- 
cept for the fact that the real forces are replaced by ef- 
fective forces defined by the direct correlation function). 

The connection between curvature and force is intu- 
itively obvious since motion on a flat manifold requires 



no force. Indeed, the connection was made by one of the 
pioneers in connecting mechanics and differential geom- 
etry, Lipschitz, who wrote in 1873 [132| : "When a mate- 
rial particle, which is not influenced by any accelerating 
force, is bound to move on a given surface, the pressure 
exerted in each point of the trajectory is inversely pro- 
portional to the r adius of curvature of this trajectory" 
(quoted from Ref. Il33l ). What happens is the following. 
When a particle (i.e., the system) moves on a perfectly 
smooth surface, since no work is performed on the parti- 
cle, the kinetic energy is conserved and thus the particle's 
velocity v is constant. This implies that the centripetal 
force keeping the particle on the surface, oc w 2 /r, is in- 
versely proportional to the local curvature radius r. 

For highly viscous liquids most motion is vibrational 
and one can estimate F by adopting a harmonic approxi- 
mation. Writing for the force on a particle — Cx where x 
is the vibrational displacement in an axis direction, one 
has (F 2 ) oc C 2 (x 2 }. Since C(x 2 )/2 = k B T/2 by equipar- 
tition, this means that if a = p -1 / 3 is the interatomic 
spacing, (F 2 ) = p- 2 / 3 (F 2 ) /[3N(k B T) 2 } cx a 2 /(x 2 ) = 
l/(x 2 ). Thus the quantity a 2 /(x 2 ) identifies tl and, in 
this sense, controls the relaxation time for s imple liq - 
uids just as entropy does. This is an old idea plll4lj . 
which is closely related to the reasoning behind the shov- 
ing model and other elastic models for the te mperatur e 
dependence of viscous liquids' relaxation time [104 \142\ - 



ll44j . Note also that, since a 2 /(x 2 ) = l/{i 2 ) identifies fi, 
for simple liquids the crossover to activated transitions 
takes place at a quasiuniversal value of the reduced-unit 
vibrational displacement (a: 2 ) as re cently predicted by 
Lubchenko and coworkers |l45l Il46l | . 

A single microconfiguration is enough to identify f2. 

0(A) quasiuniversality implies that the pair poten- 
tial plays no role for several properties of a simple liq- 
uid, namely the isomorph invariants. This is because 
given an equilibrium configuration of a simple liquid 
R = (ri, rjy), the corresponding reduced-coordinate 
vector R identifies the relevant manifold (l. From this 
the isomorph invariants like entropy, specific heat, re- 
duced diffusion constant, reduced incoherent scattering 
function, etc, can be calculated. The reduced-unit radial 
distribution function g(f) is also included in the list, of 
course, but this quantity is trivially given by R. This il- 
lustrates again the Young- Andersen approximate scaling 
principle that for simple liquids knowing g(f) determines 
many other quantities. 



IV. TOWARDS A THEORY OF SIMPLE 
LIQUIDS 

Rosenfeld in 1977 justified IPL quasiuniversality by ar- 
guing that each n-IPL system is well represented by some 
HS reference system [7J]. Since the HS system has just 
one parameter, this certainly implies IPL quasiuniversal- 
ity. However, given the fact that the HS system is the 
n — > oo limit of n-IPL systems, this reasoning may be 
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regarded as circular by assuming what is to be arrived 
at. Section [XT] gave one argument for quasiuniversality, 
but we wish here to supplement it by another argument 
suggesting an underlying reason for IPL quasiuniversal- 
ity. 

We take as starting point that the pair potential de- 
fined by a simple exponentially decaying function of r is 
strongly correlating, i.e., defines a simple liquid. This 
remains to be thoroughly investigated and docu ment ed, 
but the recent simulations by Veldhorst et al. |l47j of 
the Buckingham potential, which has a harsh exponen- 
tially repulsive term, certainly indicates that this is the 
case. Thus we assume that systems with pair poten- 
tials of the form v(r) — eexp(— r/a) are strongly corre- 
lating (at least in a large part of the (c 3 /0, e/ksT) pa- 
rameter space). This means that the reduced-coordinate 
constant-potential-energy hypersurfaces Q, which are a 
priori parameterized by the two dimensionless numbers 
a 3 p and e/fc^T, constitute a single-parameter family 
fi(A), where A = X(a 3 p,e/k B T). 

We proceed to argue that this is the quasiuniversal one- 
parameter family characterizing all simple liquids. It is 
enough to show that all IPL pair potentials have in com- 
mon these reduced-coordinate constant-potential-energy 
hypersurfaces. Now, following the reasoning of Sec. IIII1 
there is additivity among exponentially decaying pair po- 
tentials in the sense that any linear combination of two 
such potentials, C\ exp(— r/<j\) + C2 exp(— r/os), has the 
same family of reduced-coordinate constant-potential- 
energy hypersurfaces as a single exponential. By gen- 
eralization, this implies via the mathematical identity 

r- n = - / x n - x eT xr dx (12) 

0-1)! J 

that the family 0(A) is common to the IPL potentials 



and, by implication, to all simple liquids. 



V. SUMMARY 

The quasiuniversalities observed for simple liquids' 
structure, dynamics, and thermodynamics find a natural 
explanation in terms of quasiuniversality of the single- 
parameter family il(A) of reduced-coordinate constant- 
potential-energy hypersurfaces. f2(A) quasiuniversality 
was derived for IPL systems in Sec. HH but can ar- 
rived at assuming that the exponentially repulsive pair 
potential is simple (Sec. II VI) . From IPL quasiuniversality 
fi(A) quasiuniversality generalizes to all simple liquids by 
virtue of their property of having strong correlations be- 
tween NVT equilibrium virial and potential-energy fluc- 
tuations. 

It is important to emphasize that there is no exact uni- 
versality among simple liquids, only approximate univer- 
sality. A clear demonstration of this is provided by the 
well-known fact that the crystalline state is only FCC 
for IPL exponents larger than seven (below it is BCC). 
Another point to be emphasized is that, in contrast to 
mode-coupling theory and other fully renormalized the- 
ories of liquid dynamics, the present approach does not 
distinguish between short and long time scales - quasiu- 
niversality applies on all time scales (in reduced units). 



ACKNOWLEDGMENTS 

The centre for viscous liquid dynamics "Glass and 
Time" is sponsored by the Danish National Research 
Foundation (DNRF). 



[1] K. S. Pitzer, "The volumetric and thermodynamic prop- 
erties of fluids. I. Theoretical basis and virial coeffi- 
cients," J. Am. Chem. Soc. 77, 3427-3433 (1955). 

[2] P. G. Debenedetti, "Structure, dynamics and thermody- 
namics in complex systems: Theoretical challenges and 
opportunities," AICHE J. 51, 2391-2395 (2005). 

[3] B. Kirchner, "Theory of complicated liquids: Investiga- 
tion of liquids, solvents and solvent effects with modern 
theoretical methods," Phys. Rep. 440, 1-111 (2007). 

[4] J. F. Douglas, J. Dudowicz, and K. F. Freed, "Lat- 
tice model of equilibrium polymerization. VI. Measures 
of fluid "complexity" and search for generalized cor- 
responding states," J. Chem. Phys. 127, 224901-01 - 
224901-19 (2007). 

[5] B. Bagchi and C. Chakravarty, "Interplay between mul- 
tiple length and time scales in complex chemical sys- 
tems," J. Chem. Sci. 122, 459-470 (2010). 

[6] J. L. Barrat and J. P. Hansen, Basic concepts for sim- 
ple and complex liquids (Cambridge University Press, 



2003). 

[7] J. D. Bernal, "The Bakerian lecture, 1962. The structure 
of liquids," Proc. R. Soc. London Ser. A 280, 299-322 
(1964). 

[8] I. Z. Fisher, Statistical theory of liquids (University of 

Chicago, Chicago, 1964). 
[9] S. A. Rice and P. Gray, The statistical mechanics of 
simple liquids (Interscience, New York, 1965). 

[10] H. N. V. Temperley, J. S. Rowlinson, and G. S. Rush- 
brooke, Physics of simple liquids (Wiley, 1968). 

[11] Y. Rosen, "A new theory of simple classical fluids," 
Phys. Lett. A 73, 31-34 (1979). 

[12] N. K. Ailawadi, "Equilibrium theories of classical flu- 
ids," Phys. Rep. 57, 241-306 (1980). 

[13] J. S. Rowlinson and B. Widom, Molecular theory of cap- 
illarity (Clarendon, Oxford, 1982). 

[14] C. G. Gray and K. E. Gubbins, Theory of molecular 
fluids (Oxford University Press, 1984). 



8 



[15] D. Chandler, Introduction to modern statistical mechan- 
ics (Oxford University Press, 1987). 

[16] J. P. Hansen and J. R. McDonald, Theory of simple 
liquids, 3rd ed. (Academic, New York, 2005) . 

[17] F. H. Stillinger, "Phase transitions in the Gaussian core 
system," J. Chem. Phys. 65, 3968-3974 (1976). 

[18] F. H. Stillinger and D. K. Stillinger, "Negative thermal 
expansion in the Gaussian core model," Physica A 244, 
358-369 (1997). 

[19] S. Prestipino, F. Saja, and P. V. Giaquinta, "Phase 
diagram of softly repulsive systems: The Gaussian 
and inverse-power-law potentials," J. Chem. Phys. 123, 
144110-1 - 144110-9 (2005). 

[20] C. E. Zachary, F. H. Stillinger, and S. Torquato, "Gaus- 
sian core model phase diagram and pair correlations 
in high euclidean dimensions," J. Chem. Phys. 128, 
224505-1 - 224505-16 (2008). 

[21] M. Engel and H. R. Trebin, "Self-assembly of 
monatomic complex crystals and quasicrystals with a 
double-well interaction potential," Phys. Rev. Lett. 98, 
225505-1 - 225505-4 (2007). 

[22] E. A. Jagla, "Core-softened potentials and the anoma- 
lous properties of water," J. Chem. Phys. Ill, 8980- 
8986 (1999). 

[23] Z. Yan, S. V. Buldyrev, N. Giovambattista, P. G. 
Debenedetti, and H. E. Stanley, "Family of tunable 
spherically symmetric potentials that span the range 
from hard spheres to waterlike behavior," Phys. Rev. 
E 73, 051204-1 - 051204-8 (2006). 

[24] E. Lomba, N. G. Almarza, C. Martin, and C. McBride, 
"Phase behavior of attractive and repulsive ramp fluids: 
Integral equation and computer simulation studies," J. 
Chem. Phys. 126, 244510 -1 - 244510 -12 (2007). 

[25] G. Franzesc, "Differences between discontinuous and 
continuous soft-core attractive potentials: The appear- 
ance of density anomaly," J. Mol. Liq. 136, 267-273 
(2007). 

[26] A. B. de Oliveira, P. A. Netz, T. Colla, and M. C. 
Barbosa, "Structural anomalies for a three dimensional 
isotropic core-softened potential," J. Chem. Phys. 125, 
124503-1 - 124503-7 (2006). 

[27] N. M. Barraz, E. Salcedo, and M. C. Barbosa, "Ther- 
modynamic, dynamic, and structural anomalies for 
shoulderlike potentials," J. Chem. Phys. 131, 094504- 
01 - 094504-07 (2009). 

[28] N. M. Barraz, E. Salcedo, and M. C. Barbosa, "Ther- 
modynamic, dynamic, structural, and excess entropy 
anomalies for core-softened potentials," J. Chem. Phys. 
135, 104507-1 - 104507-11 (2011). 

[29] Y. D. Fomin, E. N. Tsiok, and V. N. Ryzhov, "Inver- 
sion of sequence of diffusion and density anomalies in 
core-softened systems," J. Chem. Phys. 135, 234502-1 
- 234502-8 (2011). 

[30] T. S. Ingebrigtsen, T. B. Schr0der, and J. C. Dyre, 
"What is a simple liquid?" Phys. Rev. X 2, 011011-1 - 
011011-20 (2012). 

[31] N. P. Bailey, U. R. Pedersen, N. Gnan, T. B. Schr0der, 
and J. C. Dyre, "Pressure-energy correlations in liquids. 

I. Results from computer simulations," J. Chem. Phys. 
128, 184507 (2008). 

[32] N. P. Bailey, U. R. Pedersen, N. Gnan, T. B. Schr0der, 
and J. C. Dyre, "Pressure-energy correlations in liquids. 

II. Analysis and consequences," J. Chem. Phys. 129, 
184508 (2008). 



[33] T. B. Schr0der, N. P. Bailey, U. R. Pedersen, N. Gnan, 
and J. C. Dyre, "Pressure-energy correlations in liquids. 
III. Statistical mechanics and thermodynamics of liq- 
uids with hidden scale invariance," J. Chem. Phys. 131, 
234503 (2009). 

[34] N. Gnan, T. B. Schr0der, U. R. Pedersen, N. P. Bai- 
ley, and J. C. Dyre, "Pressure-energy correlations in liq- 
uids. IV. Isomorphs in liquid phase diagrams," J. Chem. 
Phys. 131, 234504 (2009). 

[35] T. B. Schr0der, N. Gnan, U. R. Pedersen, N. P. Bailey, 
and J. C. Dyre, "Pressure-energy correlations in liquids. 
V. Isomorphs in generalized Lennard- Jones systems," J. 
Chem. Phys. 134, 164505 (2011). 

[36] T. S. Ingebrigtsen, S. Toxvaerd, O. J. Heilmann, T. B. 
Schr0der, and J. C. Dyre, "NVU dynamics. I. Geodesic 
motion on the constant-potential-energy hypersurface," 
J. Chem. Phys. 135, 104101 (2011). 

[37] R. M. J. Cotterill and J. U. Madsen, "Localized fluid- 
ity modes and the topology of the constant-potential- 
energy hypersurfaces of lcnnard-jones matter," Phys. 
Rev. B 33, 262-268 (1986). 

[38] A. Scala, L. Angelani, R. Di Leonardo, G. Ruocco, and 
F. Sciortino, "A stroll in the energy landscape," Phil. 
Mag. B 82, 151-161 (2002). 

[39] C. Wang and R. M. Stratt, "Global perspectives on the 
energy landscapes of liquids, supercooled liquids, and 
glassy systems: The potential energy landscape ensem- 
ble," J. Chem. Phys. 127, 224503-1 -224503-10 (2007). 

[40] T. S. Ingebrigtsen, S. Toxvaerd, T. B. Schr0der, and 
J. C. Dyre, "NVU dynamics. II. Comparing to four other 
dynamics," J. Chem. Phys. 135, 104102 (2011). 

[41] T. S. Ingebrigtsen and J. C. Dyre, "NVU dynamics. III. 
Simulating molecules," in preparation (2012). 

[42] T. S. Ingebrigtsen, L. B0hling, T. B. Schr0der, and 
J. C. Dyre, "Thermodynamics of condensed matter with 
strong pressure-energy correlations," J. Chem. Phys. 
136, 061102-1 - 061102-4 (2012). 

[43] N. W. Ashcroft and J. Lekner, "Structure and resistivity 
of liquid metals," Phys. Rev. 145, 83-90 (1966). 

[44] J. -P. Hansen and L. Verlet, "Phase transitions of 
the Lennard- Jones system," Phys. Rev. 184, 151.161 
(1969). 

[45] R. Grover, "Liquid metal equation of state based on 
scaling," J. Chem. Phys. 55, 3435-3441 (1971). 

[46] W. G. Hoover, S. G. Gray, and K. W. Johnson, "Ther- 
modynamic properties of the fluid and solid phases 
for inverse power potentials," J. Chem 55, 1128-1136 
(1971). 

[47] Y. Hiwatari, H. Matsuda, T. Ogawa, N. Ogita, and 
A. Ueda, "Molecular dynamics studies on the soft-core 
model," Prog. Theor. Phys. 52, 1105-1123 (1974). 

[48] S. M. Stishov, "The thermodynamics of melting of sim- 
ple substances," Sov. Phys. Usp. 17, 625-643 (1975). 

[49] Y. Rosenfeld, "Universality of melting and freezing in- 
dicators and additivity of melting curves," Mol. Phys. 
32, 963-977 (1976). 

[50] H. S. Kang, S. C. Lee, T. Ree, and F. H. Ree, "A 
perturbation theory of classical equilibrium fluids," J. 
Chem. Phys. 82, 414-423 (1985). 

[51] A. V. Indrani and S. Ramaswamy, "Universal self- 
diffusion and subdiffusion in colloids at freezing," Phys. 
Rev. Lett 73, 360-363 (1994). 

[52] A. C. Branka and D. M. Heyes, "Thermodynamic prop- 
erties of inverse power fluids," Phys. Rev. E 74, 031202- 



9 



1 - 031202-11 (2006). 

[53] D. M. Heyes and A. C. Branka, "Physical properties of 
soft repulsive particle fluids," Phys. Chem. Chem. Phys. 
9, 5570-5575 (2007). 

[54] D. M. Heyes and A. C. Branka, "Self-diffusion coeffi- 
cients and shear viscosity of inverse power fluids: from 
hard- to soft-spheres," Phys. Chem. Chem. Phys. 10, 
4036-4044 (2008). 

[55] F. de Guevara-Rodriguez and M. Medina-Noyola, "Dy- 
namic equivalence between soft- and hard-core Brow- 
nian fluids," Phys. Rev. E 68, 011405-1 - 011405-8 
(2003). 

[56] L. Lopez-Flores, P. Mendoza-Mendez, Luis E. 
Sanchez-Diaz, G. Perez-Angel, M. Chavez-Paez, 
A. Viscarra-Rendion, and M. Medina-Noyola, "Dy- 
namic equivalence between atomic and colloidal 
liquids," arXiv:1106.2475 (2011). 

[57] T. Young and H. C. Andersen, "A scaling principle for 
the dynamics of density fluctuations in atomic liquids," 
J. Chem. Phys. 118, 3447-3450 (2003). 

[58] T. Young and H. C. Andersen, "Tests of an approxi- 
mate scaling principle for dynamics of classical fluids," 
J. Phys. Chem. B 109, 2985-2994 (2005). 

[59] M. J. Pond, J. R. Errington, and T. M. Truskett, 
"Mapping between long-time molecular and Brownian 
dynamics," Soft Matter 7, 9859-9862 (2011). 

[60] M. Schmiedeberg, T. K. Haxton, S. R. Nagel, and 
A. J. Liu, "Mapping the glassy dynamics of soft spheres 
onto hard-sphere behavior," EPL 96, 36010-1 - 36010-5 
(2011). 

[61] M. G. Noro and D. Frenkel, "Extended corresponding- 
states behavior for particles with variable range attrac- 
tions," J. Chem. Phys. 113, 2941-2944 (2000). 

[62] R. J. Speedy, "The energy-density landscape for soft 
spheres," J. Phys.: Condens. Matter 15, S1243-S1251 
(2003). 

[63] H. W. Sheng and E. Ma, "Atomic packing of the in- 
herent structure of simple liquids," Phys. Rev. E 69, 
062202-1 - 062202-4 (2004). 

[64] D. M. Heyes and A. C. Branka, "The influence of poten- 
tial softness on the transport coefficients of simple flu- 
ids," J. Chem. Phys. 122, 234504-1 - 234504-9 (2005). 

[65] T. Scopigno, R. Di Leonardo, L. Comez, A. Q. R. Baron, 
D. Fioretto, and G. Ruocco, "Hard-sphere-like dynam- 
ics in a non-hard-sphere liquid," Phys. Rev. Lett. 94, 
155301-1 - 155301-4 (2005). 

[66] P. Orea, Y. Reyes-Mercado, and Y. Duda, "Some uni- 
versal trends of the Mie(n,m) fluid thermodynamics," 
Phys. Lett. A 372, 7024-7027 (2008). 

[67] D. M. Heyes and A. C. Branka, "Density and pressure 
dependence of the equation of state and transport coef- 
ficients of soft-sphere fluids," Mol. Phys. 107, 309-319 
(2009). 

[68] E. Lange, J. B. Caballero, A. M. Puertas, and M. Fuchs, 
"Comparison of structure and transport properties of 
concentrated hard and soft sphere fluids," J. Chem. 
Phys. 130, 174903-1 - 174903-8 (2009). 

[69] P. E. Ramirez-Gonzalez and M. Medina-Noyola, "Glass 
transition in soft-sphere dispersions," J. Phys.: Con- 
dens. Matter 21, 075101-1 - 075101-13 (2009). 

[70] P. E. Ramirez- Gonzalez, L. Lopez-Flores, H. Acuna- 
Campa, and M. Medina-Noyola, "Density-temperature- 
softness scaling of the dynamics of glass-forming soft- 
sphere liquids," Phys. Rev. Lett. 107, 155701-1 - 



155801-5 (2011). 

[71] T. M. Truskett, S. Torquato, and P. G. Debenedetti, 
"Towards a quantification of disorder in materials: Dis- 
tinguishing equilibrium and glassy sphere packings," 
Phys. Rev. E 62, 993-1001 (2000). 

[72] J. R. Errington, P. G. Debenedetti, and S. Torquato, 
"Quantification of order in the Lennard- Jones system," 
J. Chem. Phys. 118, 2256 (2003). 

[73] S. N. Chakraborty and C. Chakravarty, "Entropy, lo- 
cal order, and the freezing transition in Morse liquids," 
Phys. Rev. E 76, 011201-1 - 011201-9 (2007). 

[74] Y. Rosenfeld, "Relation between the transport coeffi- 
cients and the internal entropy of simple systems," Phys. 
Rev. A 15, 2545-2549 (1977). 

[75] Y. Rosenfeld, "A quasi-universal scaling law for atomic 
transport in simple fluids," J. Phys.: Condens. Matter 
11, 5415-5427 (1999). 

[76] M. J. Pond, J. R. Errington, and T. M. Truskett, "Gen- 
eralizing Rosenfeld's excess-entropy scaling to predict 
long-time diffusivity in dense fluids of Brownian par- 
ticles: From hard to ultrasoft interactions," J. Chem. 
Phys. 134, 081101-1 - 081101-4 (2011). 

[77] M. Singh, M. Agarwal, D. Dhabal, and C. Chakravarty, 
"Structural correlations and cooperative dynamics in 
supercooled liquids," J. Chem. Phys. 137, 024508-1 - 
024508-10 (2012). 

[78] R. Grover and. W. G. Hoover and B. Moran, "Corre- 
sponding states for thermal conductivities via nonequi- 
librium molecular dynamics," J. Chem. Phys. 83, 1255- 
1259 (1985). 

[79] J. J. Gilvarry, "The Lindemann and Griineisen laws," 
Phys. Rev. 102, 308-316 (1956). 

[80] A. R. Ubbelohde, Melting and crystal structure (Claren- 
don (London), 1965). 

[81] M. Ross, "Generalized Lindemann melting law," Phys. 
Rev. 184, 233-242 (1969). 

[82] F. Saija, S. Prestipino, and P. V. Giaquinta, "Evalua- 
tion of phenomenological one-phase criteria for the melt- 
ing and freezing of softly repulsive particles," J. Chem. 
Phys. 124, 244504-1 - 244504-5 (2006). 

[83] G. Malescio, P. V. Giaquinta, and Y. Rosenfeld, "Struc- 
tural stability of simple classical fluids: Universal prop- 
erties of the Lyapunov-exponent measure," Phys. Rev. 
E 61, 4090-4094 (2000). 

[84] S. Toxvaerd, unpublished. 

[85] D. C. Wallace, Statistical Physics of Crystals and Liq- 
uids (World Scientific, Singapore, 2002). 
[86] D. Bolmatov, V. V. Brazhkin, and K. Tracenko, "The 

phonon theory of liquid thermodynamics," Sci. Rep. 2, 

421-1 - 421-6 (2012). 
[87] E. N. D. Andrade, "A theory of the viscosity of liquids 

- Part 1," Phil. Mag. 17, 497-511 (1934). 
[88] G. Kaptay, "A unified equation for the viscosity of pure 

liquid metals," Z. Metallkd. 96, 1-8 (2005). 
[89] H. J. Raveche, R. D. Mountain, and W. B. Streett, 

"Freezing and melting properties of the Lennard- Jones 

system," J. Chem. Phys. 61, 1970-1984 (1974). 
[90] F. Saija, S. Prestipino, and P. V. Giaquinta, "Scaling 

of local density correlations in a fluid close to freezing," 

J. Chem. Phys. 115, 7586-7591 (2001). 
[91] J. L. Tallon, "The entropy change on melting of simple 

substances," Phys. Lett. A 76, 139-142 (1980). 
[92] P. A. Egelstaff and B. Widom, "Liquid surface tension 

near the triple point," J. Chem. Phys. 53, 2667-2669 



10 



(1970). 

Y. Rosenfeld, "New method for equation-of-state cal- 
culations: Linear combinationsws of basis potentials," 
Phys. Rev. A 26, 3633-3645 (1982). 
S. A. Khrapak and G. E. Morrill, "Accurate freezing 
and melting equations for the Lennard- Jones system," 
J. Chem. Phys. 134, 094108 (2011). 
C. Alba-Simionesco, D. Kivelson, and G. Tarjus, "Tem- 
perature, density, and pressure dependence of relax- 
ation times in supercooled liquids," J. Chem. Phys. 116, 
5033-5038 (2002). 

F. Sciortino, W. Kob, and P. Tartaglia, "In- 
herent structure entropy of supercooled liquids," 
|Phys. Rev. Lett. 83, 3214-3217 (1999)] 

B. Doliwa and A. Heuer, "Finite-size effects in a su- 
percooled liquid," J. Phys.: Condens. Matter 15, S849- 
S858 (2003). 

Q. Yan, T. S. Jain, and J.J. de Pablo, "Density- 
of-states Monte Carlo simulation of a binary glass," 
|Phys. Rev. Lett. 92, 235701-1 - 237501-4 (2004)"] 
Y. Gebremichael, M. Vogel, M. N. J. Bergroth, F. W. 
Starr, and S. C. Glotzer, "Spatially heterogeneous dy- 
namics and the Adam-Gibbs relation in the Dzugutov 
liquid," J. Phys. Chem. B 109, 15068-15079 (2005). 
U. R. Pedersen, T. B. Schr0der, and J. C. Dyre, "Repul- 
sive reference potential reproducing the dynamics of a 
liquid with attractions," Phys. Rev. Lett. 105, 157801-1 
-157801-4 (2010). 

T. S. Ingebrigtsen and J. C. Dyre, unpublished. 

C. A. Angell, "Formation of glasses from liquids and 
biopolymers," Science 267, 1924-1935 (1995). 

G. P. Johari, "Introduction to the glassy state in the 
undergraduate curriculum," J. Chem. Educ. 51, 23-27 
(1974). 

J. C. Dyre, "The glass transition and elastic models 
of glass-forming liquids," Rev. Mod. Phys. 78, 953-972 
(2006). 

D. Kivelson, G. Tarjus, X. Zhao, and S. A. Kivelson, 
"Fitting of viscosity: Distinguishing the temperature 
dependences predicted by various models of supercooled 
liquids," Phys. Rev. E 53, 751-758 (1996). 

C. De Michele, F. Sciortino, and A Coniglio, "Scaling in 
soft spheres: Fragility invariance on the repulsive poten- 
tial softness," J. Phys.: Condens. Matter 16, L489-L494 
(2004). 

K. Niss, C. Dalle-Ferrier, G. Tarjus, and C. Alba- 
Simionesco, "On the correlation between fragility and 
stretching in glass-forming liquids," J. Phys.: Condens. 
Matter 19, 076102-1 - 076102-25 (2007). 
R. W. Zwanzig, "High-temperature equation of state by 
a perturbation method. I. Nonpolar gases," J. Chem. 
Phys. 22, 1420-1426 (1954). 

H. C. Longuet-Higgins and B. Widom, "A rigid sphere 
model for melting of argon," Mol. Phys. 8, 549-556 
(1964). 

B. Widom, "Intermolecular forces and the nature of the 
liquid state," Science 157, 375 (1967). 

D. Weeks, D. Chandler, and H. C. Andersen, "Role of 
repulsive forces in determining the equilibrium structure 
of simple liquids," J. Chem. Phys. 54, 5237-5247 (1971). 
S. Torquato and F. H. Stillinger, "Jammed hard-particle 
packings: From Kepler to Bernal and beyond," Rev. 
Mod. Phys. 82, 2633-2672 (2010). 



[113; 

[114 

[115; 
[ii6; 

[117 

[us; 

[119; 
[120; 

[121 

[122; 
[123; 
[124; 

[125; 

[126 
[127 

[128; 
[129; 

[130; 

[131 

[132; 
[133; 
[134; 



C. Brito and M. Wyart, "Geometric interpretation of 
previtrification in hard sphere liquids," J. Chem. Phys. 
131, 024504-1 - 024504-14 (2009). 
Y. Rosenfeld, "Self-consistent density functional theory 
and the equation of state for simple fluids," Mol. Phys. 
94, 929-935 (1998). 

A. B. Bestul and S. S. Chang, "Excess entropy at glass 
transformation," J. Chem. Phys. 40, 3731-3733 (1964). 

G. Adam and J. H. Gibbs, "On temperature dependence 
of cooperative relaxation properties in glass-forming liq- 
uids," J. Chem. Phys. 43, 139-146 (1965). 

J. Dudowicz, K. F. Freed, and J. F. Douglas, "General- 
ized entropy theory of polymer glass formation," Adv. 
Chem. Phys. 138, 125-222 (2008). 

J. C. Dyre, T. Heckhsher, and K. Niss, "A brief cri- 
tique of the Adam-Gibbs entropy model," J. Non-Cryst. 
Solids 355, 624-627 (2009). 

M. Dzugutov, "A universal scaling law for atomic diffu- 
sion in condensed matter," Nature 381, 139-141 (1996). 
N. J. Hicks, Notes on Differential Geometry (van Nos- 
trand Reinhold, New York, 1963). 

P. Dombrowski, "Krummungsgrossen gleichungs- 
definierter Untermannigfaltigkeiten Riemannscher 
Mannigfaltigkeiten," Math. Nachr. 38, 133-180 (1968). 
S. Gallot, D. Hulin, and J. Lafontaine, Riemannian 
Geometry (Springer, Berlin, 2004). 

M. Berger, A Panoramic View of Riemannian Geometry 
(Springer, Berlin, 2003). 

R. Goldman, "Curvature formulas for implicit curves 
and surfaces," Comput. Aided Geom. Des. 22, 632-658 
(2005). 

Y. Giga, Surface evolution equations: A level set ap- 
proach (Birkhauser, Basel, 2006). 

J. Bigeleisen, M. W. Lee, and F. Mandel, "Mean square 
force in simple liquids and solids from isotope effect 
studies," Acc. Chem. Res. 8, 179-184 (1975). 

H. H. Rugh, "Dynamical approach to temperature," 
Phys. Rev. Lett. 78, 772-774 (1997). 

J. G. Powles, G. Rickayzen, and D. M. Heyes, "Tem- 
peratures: old, new and middle aged," Mol. Phys. 103, 
1361-1373 (2005). 

K. S. Schweizer, "Relationships between the single par- 
ticle barrier hopping theory and thermodynamic, dis- 
ordered media, elastic, and jamming models of glassy 
systems," J. Chem. Phys. 127, 164506-1 - 164506-9 
(2007). 

M. Tripathy and K. S. Schweizer, "The influence of 
shape on the glassy dynamics of hard nonspherical parti- 
cle fluids. II. Barriers, relaxation, fragility, kinetic vitri- 
fication, and universality," J. Chem. Phys. 130, 244907- 
1 - 244907-09 (2009). 

M. Tripathy and K. S. Schweizer, "Activated dynam- 
ics in dense fluids of attractive nonspherical parti- 
cles, ii. elasticity, barriers, relaxation, fragility, and 
self-diffusion," Phys. Rev. E 83, 041407-1 - 041407-9 
(2011). 

R. Lipschitz, "Extrait de six memoires publies dans le 
Journal de Mathematiques de Borchardt," Bull. Sci. 
Math, et Astron. 4, 297-320 (1873). 
J. Liitzen, "Interactions between mechanics and differ- 
ential geometry in the 19th century," Arch. Hist. Exact 
Sci. 49, 1-72 (1995). 

C. P. Flynn, "Atomic migration in monatomic crystals," 
Phys. Rev. 171, 682-698 (1968). 



11 



[135] R. W. Hall and P. G. Wolynes, "The aperiodic crystal 
picture and free energy barriers in glasses," J.Chem. 
Phys. 86, 2943-2948 (1987). 

[136] U. Kohler and C. Herzig, "On the correlation be- 
tween self-diffusion and the low-frequency LA 2/3(111) 
phonon mode in bcc metals," Philos. Mag. A 58, 769- 
786 (1988). 

[137] U. Buchenau and R. Zorn, "A relation between fast and 

slow motions in glassy and liquid selenium," Europhys. 

Lett. 18, 523-528 (1992). 
[138] A. Heuer and H. W. Spiess, "Universality of the glass 

transition temperature," J. Non-Cryst. Solids 176, 294- 

298 (1994). 

[139] A. P. Sokolov, A. Kisliuk, D. Quitmann, A. Kudlik, and 
E.Rossler, "The dynamics of strong and fragile glass 
formers: Vibrational and relaxation contributions," J. 
Non-Cryst. Solids 172, 138-153 (1994). 

[140] D. S. Sanditov and S. S. Sangadiev, "Condition of glass 
transition in the fluctuation free volume theory and the 
Lindemann criterion for melting," Glass Phys. Chem. 
24, 285-294 (1998). 

[141] F. W. Starr, S. Sastry, J. F. Douglas, and S. C. 
Glotzer, "What do we learn from the local geometry of 



glass-forming liquids?" Phys. Rev. Lett. 89, 125501-1 - 
125501-4 (2002). 

[142] J. C. Dyre, T. Christensen, and N. B. Olsen, "Local 
elastic expansion model for viscous-flow activation en- 
ergies of glass-forming molecular liquids," Phys. Rev. B 
53, 2171-2174 (1996). 

[143] L. Larini, A. Ottochian, C. De Michele, and D. Lep- 
orini, "Universal scaling between structural relaxation 
and vibrational dynamics in glass-forming liquids and 
polymers," Nat. Phys. 4, 42-45 (2008). 

[144] J. C. Dyre and W. H. Wang, "The instantaneous shear 
modulus in the shoving model," J. Chem. Phys. 136, 
224108-1 - 224108-5 (2012). 

[145] V. Lubchenko and P. G. Wolynes, "Theory of struc- 
tural glasses and supercooled liquids," Annu. Rev. Phys. 
Chem. 58, 235-266 (2007). 

[146] P. Rabochiy and V. Lubchenko, "Liquid state elasticity 
and the onset of activated transport in glass formers," 
J. Phys. Chem. B 116, 5729-5737 (2012). 

[147] A. A. Veldhorst, L. B0hling, J. C. Dyre, and T. B. 
Schr0der, "Isomorphs in the phase diagram of a model 
liquid without inverse power law repulsion," Eur. Phys. 
J. B 85, 21 (2012). 



